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Abstract 

We describe the distribution of the first finite number of eigenvalues in a newly-forming band of the spectrum 
of the random Hermitean matrix model. The method is rigorously based on the Riemann-Hilbert analysis of 
the corresponding orthogonal polynomials. We provide an analysis with an error term of order N^'^^ where 
1/7 = + 2 is the exponent of non-regularity of the effective potential, thus improving even in the usual case the 
analysis of the pertinent literature. 

The behavior of the first finite number of zeroes (eigenvalues) appearing in the new band is analyzed and 
connected with the location of the zeroes of certain Freud polynomials. In general all these newborn zeroes 
approach the point of nonregularity at the rate N whereas one (a stray zero) lags behind at a slower rate of 
approach. The kernels for the correlator functions in the scaling coordinate near the emerging band are provided 
together with the subleading term: in particular the transition between K and K + I eigenvalues is analyzed in 
detail. 
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1 Introduction 

In this paper we consider the Hermitean matrix model in the scaling regime or -which is the same- the orthogonal 
polynomials on the real line with a varying weight, in the same spirit as [1, 2]. We address a particular situation 
of "nonregular" or "critical" potential: this means that the mean-field electrostatic potential vanishes at some 
point ^0 outside of the support of the equilibrium measure. This situation corresponds to a recent investigation 
[3] and is the situation where a band in the spectrum of the corresponding matrix model or a new component of 
accumulation of the zeroes of the orthogonal polynomials is about to appear (or has just disappeared) . 

In the paper we will be mostly taking the point of view of approximation theory and hence we will focus on 
the orthogonal polynomial side, but the more physically-oriented reader will have no difficulty in translating those 
results; in a picturesque way we will think of the zeroes of the polynomials as a growing population and thus call 
"colonies" the first zeroes appearing near the new band, also termed the "outpost" . This should explain the catchy 
title. 

We will be using a particular (simplified) version of the double-scaling limit: in this approach we keep the 
potential V{x) and the total charge T fixed but we add a piecewise constant perturbation of order ln(7V) /N to 
the potential near the outpost. In due time (App. A) we will explain how this simplified approach yields in fact 
identical results to the usual double-scaling. 

In order to explain in more detail the setup, suppose that the effective potential [4] (p{x) — ^{V{x) — 2g + £) 
vanish at as C{x — ■^o)^''^^ (with some C > 0); here is some point outside of the support of the equilibrium 
measure [4]. We then modify the potential (Sect. 3) by adding a step- like perturbation of the form 

Vix) ^ Vix) = V{x) - In NxAx) , 7 := ■ (1-1) 

Here xj is the characteristic function of a small interval around ^o- The real parameter x determines the strength 
of the perturbation and the constants are crafted for later convenience. 
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For X < the orthogonal polynomials do not exhibit any peculiar behavior in the large limit; for positive 
values of >c new zeroes of the OPs start appearing near the outpost. It is natural that -since we can have 
only an integer number of such roots- there are transition points in the asymptotic behavior for special values 
of XT. Specifically the normalizations we have chosen are such that in the asymptotic regime there are K roots 
near the outpost, where K is the integer nearest to xr. Clearly transitions must occur at >c € N + ^. The 
phenomenon is already captured by the leading order asymptotics (Sect. 3): indeed one can construct a uniform 
approximation to the orthogonal polynomials to order N~'^^^~'^^^^\ The approximation is best when >c is an integer 
and gets progressively worse and worse as >f approaches a half-integer. In particular when >r is a half integer the 
approximation breaks down (the error term is no longer vanishing as N — s- oo). This apparent obstruction was 
noted in [3] (see also [5, 6]). 

We then show how to obtain an improved approximation in Sect. 4; indeed we construct a uniform asymptotic 
solution with an error term of order N~'^'' (uniformly in x!). Even in the ordinary case x = (which was dealt 
with in [2]) our approximation is better than the one usually provided in the literature (in fact in Sect 4.7 we show 
how to obtain an approximation of order N^^ for arbitrary >r). 

Using this information we can study in detail the asymptotic behavior of the first roots (Sect. 4.4); in a 
scaling parameter C, ~ N'^ {z — ^q) they are related to the roots of the Freud's polynomials for the weight e~^ ^ dC; 
in particular 

• if X e [K,K + 1/2) then they are within a distance 0{N~''^''^) from the roots (in the C^coordinate) of the 
K-th Freud polynomial. 

• if X e (K — 1/2, K) then there are still K roots; however while K —1 of them are within 0{N~^^^) from the 
roots of the {K — l)-st Freud polynomial, the "last one" lags behind and meanders at a distance 0{N'^'^^) 
in the C-coordinate. Note that -while escaping to infinity in the ^-scaling parameter, such root is actually 
converging to the outpost at a rate N~'^^'^'^^ . 

• when X G N + 1/2 (namely (5 = ^ ) then there is a stray root that remains at a finite distance from the 
outpost, while the remaining converge to it. 

In the main body of the paper we make the simplifying assumption that the support of the equilibrium measure 
consists of one interval (one— cut assumption). However this is only a simplification and none (or almost) of the 
conclusions are at all dependent on it, but of course the formulae for the outer parametrix are much simpler to 
write and easier to handle also for those readers who do not know well the theory of Theta functions. 

In App. B we show (in a somewhat sketchy form) how to generalize to an arbitrary number of cuts: only one 
detail cannot be fully addressed in this general case, and concerns with the improved asymptotic for some of the 
exceptional values of x G N + 1/2 and exceptional spectral curves, namely the equivalent of formula 4.51 for the 
multi-cut case and the corresponding sign. The knowledge of the sign of (the imaginary part of) (4.51) is necessary 
to ensure that in particularly exceptional circumstances certain denominators (4.22) do not vanish. 

Remark 1.1 In a strange twist of events while the present manuscript was in the latest phases oj preparation, a 
similar preprint [5] has appeared where the author uses a RH analysis for the simplest nonregular case. Shortly 
(two days) after another independent preprint [6] on the same topic has appeared, dealing with a more general type 
of nonregularity (of the same type we deal). 
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Our work however provides a refined error analysis up to order whereas both the previous papers apparently 

give only the leading term asymptotics, with an error term of order _/V'^')'(i+2'5)^ thus not valid at the transition points 
(although in [5] an analysis of the half-integer case is also provided but only for the simplest nonregularity exponent) . 

In addition we provide a detailed analysis of the location of the zeroes of the orthogonal polynomials near the 
outpost. 



2 General setting 



The setting of the present paper will be identical in the most part to [2], [1]. Although we present in a self-contained 
way we refer the reader to the pertinent literature for more details on the subject. 
Consider the Hermitean matrix model with measure given by 

1 



(2.1) 



where Zj^ is a normalization constant. It is known [7] that the model is "solvable" in terms of orthogonal polyno- 
mials (OPs) and that all spectral statistics can be described in terms of suitable kernels constructed in terms of 
OPs. 

Let {pn{x)} be the corresponding (monic) OPs that satisfy the following orthogonality condition with the 
potential V{x) which we assume to be real and analytic. 



The spectral statistics of the model is determined by the ChristofFel-Darboux kernel [7] 



3=0 



Pj{x)pj{x') _ Pn{x)pn-l{x') - Pn-l{x)pn{x') 



(2.2) 



(2.3) 



hj hn-i{x — x') 

The OP are uniquely characterized by the following Riemann-Hilbert problem. Define for z G C \ M the matrix 

p„(x)e^T^(^) dx 



Y{z) := y„(z) := 



Pn{z) 



0n(^) 



<Pn{z) 



1 

2iTT 



(2.4) 



The above matrix has the following jump-relations and asymptotic behavior that uniquely characterize it 
[8, 9, 10, 11] (we drop the explicit dependence on n for brevity) 



F+(x)-F_(x) 



1 e-T^(==) 
1 



Y{z)^{l + 0{z-')) 



z" 




(2.5) 



Replacing the orthogonality condition (2.2) by the above jump (and boundary) conditions (2.5) we obtain the 
Riemann-Hilbert problem for the OPs. Using this setup we especially want to investigate the asymptotics of the 
OPs as their degree n := N + r goes to infinity while r being fixed to an integer. 

The first step to solve the problem is to find the g-function, which we will define using the equilibrium 
measure below. We briefly recall that, in the simple case where the contour of integration in (2.2) is the real axis 
(see [12, 13] for a general approach not relying on a variational problem), the equilibrium measure is obtained from 
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the solution of a variational problem for a functional over probability measures on the real axis, in the sense of 
potential theory [4]. Indeed define the weighted electrostatic energy [4] 

J'[tj]:^2 [ V{x)dfi{x)+ [ [in ^ dfi{x) dfi{x') (2.6) 

where d/i is a positive measure supported on the real axis with total mass T ~ d/i(x). 

It is known that the functional J-' attains a unique minimum (under mild assumptions on the growth of V(x) 
at infinity) at a measure p that is called the equilibrium measure [4, 15]. 

It is also known [14] that the support of the measure p consists of a finite union of disjoint bounded intervals 
and that p is smooth on the interior of the support. 

Taking avail of the equilibrium measure, the g function [15] is then typically defined as 

g{z):^ I p{x)ln{z- x)dx^Tlnz + 0{z-^), (2.7) 

where the logarithm must be defined with an appropriate cut extending -say- from the leftmost endpoint of the 
support of p to +00. The derivative of the function g{z) (the "resolvent") satisfies a pseudo-algebraic equation 
which is key to many considerations in a different context but will be mostly irrelevant in this paper. 

The main properties that enter the steepest descent analysis are the standard properties of the logarithmic 
transform. To this end we note that the representation (2.7) implies immediately that ^g{x) is harmonic away 
from the support of p and continuous on the whole complex plane. The Euler-Lagrange variational equations 
equivalent to the optimality of the equilibrium measure p [4] can be rephrased in terms of the following conditions 
for the (^-function 

• for X CzR we have 

%(x)>0, ^{z):^^-g{z)+^- = ^- J p{y)Hx~-y)dy+^- (2.8) 

for a suitable real constant £. di(p is the effective potential of the related electrostatic problem. 

• The opposite inequality (and hence the equality) holds on the support of p. Especially, £ is chosen such 
that ^(f = on the support of p. (The support of p will be called the cuts because they form the cuts of 
the functions g'{z) and (p'{z).) Here and in the previous point, the g-function should be understood as the 
analytic function defined by its integral representation (2.7) on the simply connected domain obtained by 
removing a half-line starting e.g. at the rightmost endpoint of the support of p and extending towards — c». 
Then ^(p(x) is actually nothing but the boundary-value l/2{(p^{x) + ip^{x)) 

• In suitable finite left /right neighborhoods of the cuts, the function ^ip(x), which is also harmonic in the 
domain of analyticity of V{x), is negative. 

The situation we want to address in this paper is the case where the inequality (2.8) is not strict on M outside 
the cuts, but at some point outside of the cuts, the inequality is an equality. This situation corresponds to the 
critical situation where a cut is about to emerge (or has just disappeared) at ^o- In sai optimistic view, we will 
look at the situation as of that of an emerging spectral band, being gradually populated by eigenvalues; it is thus 
appropriate to refer to the neighborhood of the emerging band as an outpost colony of eigenvalues. 
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The pair of potentials V{x) and total charge T for which the inequality (2.8) is not strict either outside or 
inside the cuts (this last occurrence corresponding to the merger of two cuts) are called nonregular or, in the 
more physical oriented literature, critical. The steepest descent analysis was completely carried out in [2] in 
general terms and specific analysis linking with Painleve theory was carried out in [1] in the case of a merger. For 
the "birth of a cut" (i.e. outpost colonization), a heuristic arguments and a double-scaling approach'^ were used 
in [3]; while the intents of our note and of [3] are clearly the same, the methods employed are radically different. 

The conclusion that we achieve in rigorous mathematical way will be -however- parallel to that of [3] , namely 
to show that we can describe the statistics of the first finite number of eigenvalues that populate the forming cut 
in terms of an effective "microscopic" matrix model of the size of the population of the outpost. In addition we 
strengthen those result by localizing exactly the relevant, finite number of roots of the orthogonal polynomials. 



3 Modified setting: changing chemical potential 



In order to study a nontrivial scaling limit we modify 
the setting of the problem as follows. As mentioned in 
the previous section, we consider the situation when we 
have 5R(^(^o) = for G M outside the cuts. Suppose 
that ^0 is the point where ip{z) has the critical point of 
order 2i/ + 2, i.e. 



^(z) = o((2-eo)'''+2), 



Co. 



(3.1) 



Then we choose a finite open interval J containing 
^0 that does not contain any other turning points. We 
then consider the following modified orthogonality 
relation: 



+iV2-^ I p, 



^{x)pm{x)e "^("^Mx, 



(3.2) 



where we have defined the exponent of nonregular- 

ity 

7 (2iy + 2)-i. (3.3) 




Figure 1: The potential V with the chemical potential 
added, and the corresponding effective potential. 



In the above relation the parameter x e M will eventually determine the size of the population of the colony near 
Co. 

As the reader may realize this amounts simply to a step-wise modification of the potential: if x.J is the 
characteristic function of the interval J then we may rewrite (3.2) as a single integral without A'^^^t using the 
modified potential 



V{x) := V{x) 



lnN]xj{x). 



(3.4) 



^This means that not only n —> oo but also V, T are let depend on n in a "slow" and fine— tuned way. 



6 



The most interesting regime will turn out to be x > 0, so that the potential is slightly depressed near (Fig. 1). 
While this is a "small" perturbation of the potential (which would be irrelevant in a noncritical situation), since 
our potential is critical the effect of this perturbation is fine-tuned to obtain a nontrivial perturbation. It will also 
become clear that the actual choice of J is irrelevant as long as it contains ^ind no cuts. 

Although this "discontinuous" deformation may seem quite artificial at first, it should become apparent later on 
that it actually makes no difference on the actual behavior near the outpost. In a certain sense this is the essence 
of universality, but we will explain in Appendix A how to approach the same problem from a more "canonical" 
double-scaling limit, while retaining the main features. 

The advantage of this simplified approach is that allows us to immediately concentrate on the significant features 
(the actual RHP) without hindering the analysis into details regarding the appropriate ^-function. 

3.1 Normalized and lens-opened RHP 

Taking avail of the general wisdom, in order to streamline the derivation we open the lenses before normalizing 
the problem^, thus modifying the jumps as shown in figure 2. Lens opening simply means that we redefine 



^ew 
5^ew 



Y 
Y 



-CT 
1 



on the upper lip, 
on the lower lip. 



(3.5) 
(3.6) 



1 n 



1 




^ * 7\ 1 

















Figure 2: The jump matrices for Y. 
After the lens-opening we define, 



Y{z) e^''-'Y{z)e 
which satisfies a new, simpler RHP: 

Y{x)+ = Y{x) 
Y{z) ~ (H-0(z"i))z'''"^ 



- 5T3(2)cr3 



Q 2T 



For the time being the jumps 
on the green circles are the iden- 
tity but later we define separate 
RHP problems inside the circles. 
Then we will call all the RHP 
inside one of the green disks the 
local problem whereas we call 
the problem outside the outer 
problem. 

(3.7) 



ef (9--9+) g-f (V-g+-s_+«)^2>r7x.,(z) 



(3.8) 
(3.9) 



For simplicity we assume that there is only one finite band in the spectrum, namely the spectral curve is of 
genus 0; the generalization to more bands is not conceptually a problem but requires the use of 0~-functions which 



We are of course assuming that V{z) is real-analytic. 
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would make the note quite more technical and long. Under this assumption, for x e M, 

g+{x) ~ —g^{x) + V{x) + I , for a; e M on the cut, 
g+{x) = g-{x) — 2inT , for x e M on the right of the cut, 
g+{x) ~ g-{x) , for a; e M on the left of the cut. 



(3.10) 
(3.11) 
(3.12) 



Everywhere else on the complex plane g{z) is holomorphic. On account of these properties for the (^-function the 
jumps for Y are shown in the figure 3. 

In the following the size of 
the green circles will be fixed to 
a nonzero value. In this case the 
reader can verify that -outside 
of the green disks- the jumps on 
the black and red lines become 
exponentially close to the iden- 
tity, and uniformly so in C] 
L°°. 




Figure 3: The jump matrices for Y . 



3.2 Outer parametrix 

For simplicity of exposition we assume that there is only one spectral band (1-cut) apart from the one that is about 
to emerge. Furthermore we assume that the irregular point of the potential problem is set to = 0; without loss 
of generality. 

Removing all the jumps that 
are exponentially close to the (."i o) 

identity, we are left with the ( ] ^ ( )_ 

jump matrix as shown in Figure 
4. This provides the asymptotic 
RHP that we will use to define 
outer parametrix. Below we describe the RHP that the outer parametrix satisfies. (For a specific solution to 
the RHP we will use with a subscript or a dressing such as in ^i<-.) 




Figure 4: Jump matrix for 



*(z) ~ (1 + 0(z-i))z'^'^% 
1 



-1 



~ oo, 

on the cut. 



It also needs to be supplemented by the boundary conditions at the turning points, 

*(z) = 0{{z - a)"j), *(z) = 0{{z - b)-i), 



(3.13) 
(3.14) 

(3.15) 



where a and b are the two turning points. For the specific outer parametrix that we consider soon, the growth 
condition at the outpost is given by 



A^{z) B^{z) 
Ay{z) By{z) 



(3.16) 
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The matrix [A/f,Bif] is analytic at z (note that — 0) and det[Aif,BA'] = 1- The above four conditions give 
the Rieniann-Hilbert problem for ^P. We now describe a specific solution (and [Ax,Bx]) to the above RHP. 

We now use the fact that there is only one cut so that the two-sheeted cover of the z-plane is a rational (i.e. 
genus 0) curve. The modifications needed for the case of an arbitrary number of cuts are sketched in appendix. 

Let t be the uniformizing map of the genus-0 Riemann surface. We let on the i-plane to map to the outpost on 
z-plane. For the simplicity of the normalization we choose the location of the cut and the outpost in the following 
way. 



z{i) := 



6- 



t + 



1 



6+a a — b , s/1 



(3.17) 



4 \ t J 2 

where a < b are the endpoints of the band; in the i-plane they correspond to t — ±1. There are 2 points in the 
t-plane projecting to z = namely the two solutions of z{t) = 0. We denote the one outside the unit circle by tQ, 
the other being ^. 

Define the spinorial Baker Akhiezer vectors 



^'(i)(t) := 



-if 

ir 



t-l/to 



K 



dt 



1 / t-to \ 
1^ t-l/to J 

■-1 ( t-l/to A 



dt 



r 



t-l/to 

t-t 



dt 
di 



3.18) 




to 



i-plane 



where K is the closest nonnegative integer to x 
(e.g. if X = 2.4 then i^T = 2, if >r = 2.6 then K = 3). 
Note that the definition is ambiguous for x e ^ + Z; 
indeed it will be seen that for these exceptional val- 
ues we cannot obtain a strong asymptotic result using 
these methods and we need to use a refinement (Sect. 
4), and the asymptotic for of the OPs has a discontin- 
uous change, namely the model exhibits a nonlinear 
Stokes' phenomenon in >c. 

The advantage of this spinor representation and the 
uniformizing coordinate is that we can easily write a 
general solution to the RHP (3.13),(3.14),(3.15),(3.16) ^-^^^^ ^. r^^^ uniformization of the plane shced along the 
of the outer parametrix when there are exceptional support of the equihbrium measure. The pattern region 
points such as the outpost. is the "unphysical sheet". 

Using (3.18) one can write the following solution. 



- •- 

Co 



z-plane 



(^) 



rtT3 + i 



'dz 



vl/«(i(z)),vl/(i).*(i(z)) 



jb-a) 
4 



4^'(t) 
b—a 



t-tp 

t-l/tv. 



(3.19) 

t=t{z) 



In the above t{z) is the determination of t that lies outside the unit circle for z not on the cut (the physical sheet) 
and the appropriate value on the cut viceversa. The squareroot is the determination that behaves as 
near t = oo. 



b—a 



1 
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We see also that the above matrix behaves as follows for large x and has the following jump-discontinuity 

(3.20) 



1 



-1 



It is to be noted that 



^k{z) = 0{{z - a)-i) , ^k{z) = 0{{z - . 
More importantly, the parametrix we have constructed behaves as follows in the vicinity of the outpost 

vI/k(z) = (C + 0(z))z^^^ 

From the expression of 'i'K (3.16) one can also write [Ai^,Bx] as follows 



[Ak,Bk] = 



[Ak(0),Bk(0)] + [Ai,(0),B'^(0)]z + 0{z^), 



where we write the first two terms in the expansion for future reference. 



[A^(0),B;^(0)] 



[A^(0),B'^(0)] = 



\ 4 ) ^0 



1 - in 



it 



''0 



\ 4 ) ^0 



1 - in 



t() V ab 





b 



[Aa(0),Ba(0)] 



Aab 



(73 



KlT3 



K{b+a) 



(ig - l)Vab) Vab"'" ' 2ab 



-0-3 



(3.21) 
(3.22) 

(3.23) 
(3.24) 

(3.25) 
(3.26) 
(3.27) 



Remark 3.1 In [6, 5] the authors used a scaling limit by removing the charges corresponding to the new zeroes by 
adding a point-wise charge with the same total mass. This amounts to multiplying the jump by a factor (a;— ^o)^^ ( 
our notation). As a drawback they need to cure the non-constant jump residual after the lens-opening by introducing 
a scalar function D(x) (Szego function) solving a new (scalar) RHP on the cut. In our case this scalar function is 

"built-in" the outer parametrix and corresponds to the term {^-j^^^^ . As a result the outer parametrices are at 
first sight of different nature, but -as it should- all terms can be put in correspondence in the three approaches. 

3.3 Parametrix near the simple turning points 

This part is essentially identical to the established results in [15, 2]. We con- 
sider only the right turning point at z = 6; the method works for the other 
one. Inside the green disk we solve the exact RHP. 

For simplicity of exposition we consider the case where the turning point 
is simple, namely 



^{z) = -V{z) - g{z) = C{z - &)5(1 + 0{z - &)), 



(3.28) 



A2 

1 
-1 



1 Q-riV+g^+g- 
1 



^3 



for z not on the cut. [We will indicate the trivial modifications needed in case 
of nonregular behavior later.] 
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Figure 6: The jumps of the exact 
solution Y near the soft-edge. 



3 N . . 

2 —ip{z), 



(3.29) 



We define a local coordinate by 

2 

3^ ■ T 

where the determination of the root is such that the (blue) cut is mapped to 
]R_ of (-plane. 

We then introduce the standard Airy parametrix A'^{£,) [15, 2] as the piecewise defined matrix (see Figure 
6) constructed in terms of the Airy function Ai(x) as follows 



-4? (6 



2/0 

y'o 



-y2 
-y'2 



ea 



3 

C2 0-3 



27re < 



-2/1 
-2/i 



"^Mei" 
-2/2 



3 

C2<T3 



= 1 



= 2 



-2/2 2/i)g|cS.3 
-2/2 2/1 

2/0 2/1 



2/0 2/1 



5 3 



where we have used the definitions, 
Each of the above has the following uniform asymptotic behavior near ^ ^ oo. 

^''(c)-r-^(_\ ;)e-5-(i + o(r^/^)). 

This matrix has jumps on the rays (in the ^-plane) as indicated in Fig. 7 below. 

Thus, the final form of the local parametrix is simply 

■.=Fiz) 




e 4 



1 

71 



1 -1 
1 1 



Such matrix has the properties 

• It solves the exact jump conditions of the RHP in Fig. 6; 

• the prefactor F{z) solves a RHP on the left 



Fiz)+ = 



-1 

1 



F{z)- 



Figure 7: The usual jumps for the 
local parametrix 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



(3.34) 



• it behaves as .4(^) = 1 + 0{N^^) uniformly on the boundary on account 
that z^"'^ is analytic and invertible with analytic inverse in the neigh- 
borhood. The only point to raise is that 'i>{z)A{z) is bounded inside 
the disk; indeed near the turning point we have 



«' = C'((z-c)~i), F(z) = C'((z-c)-J). 



(3.35) 
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Thus the product ^>(z)F(z) may at most have squareroot singularities: however, comparing the RHP that 
they solve, we see that the product is a single- valued matrix, thus must be analytic since at worst it may 
have singularities of type (z — c)^5. Since F{z) is solely responsible for the singularities arising in the local 
parametrix a,t z = b, this proves the assertion. 

Remark 3.2 If the turning points are non-regular then we should use the local parametrices for the usual problem 
as described in [2]. None of the above considerations (except for the bounds in the error terms) is significantly 
modified. 



3.4 The local parametrix at the outpost 



Viz) 



At the outpost the effective potential behaves as (p{z) — ^ 
new conformal parameter z as follows: 

1 



-z-.^CP-^izy = z + 0{z') 



TCoz^''+^ with Co > 0. We define a 



(3.36) 



We define D to be a finite open neighborhood around z — Q that maps univalently to a disk centered at z — Q. We 
also define the local coordinate 

( := C2N^~z. (3.37) 

The RHP satisfied by the local parametrix R is as follows. (We will use R with a subscript or a decoration such as 
in Rk, to indicate a specific solution to the RHP.) 



1 e-^""^^ N'^'^^ 



1 

R - 1 + 0(7V^') , zee 



(3.38) 



(3.39) 



Here e is some positive number that will be determined in the subsequent analysis. Increasing e leads to a better 
asymptotics. In addition to the above conditions, we also require that ^'i? is analytic in D, which implies that 
deti?= 1. 

In section 3.2, we obtain an outer parametrix k which has the pole behavior of order K at the outpost. We 
first look for the corresponding local parametrix, which we will call Rk- 

3.4.1 Case x < 

We first observe that if >r < 0, then the solution is immediately written as^ 

Ro:= (3.40) 



1 
On i9D we have 

i? - 1 + ©(TV^-^^-^) . (3.41) 

This situation is "trivial" from the point of view of the asymptotics and hence we will only focus on the case x > Q 
in the following. 

^In fact all the results hold for < 1/2 as well. 
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Remark 3.3 In fact we observe from the explicit solution that ij x < —1v — 2 then we may simply use the 
identity matrix for the local parametrix, thus committing an error smaller than 0{N~^) which anyway arises on 
the boundary of the other turning points. 

3.4.2 Case x > 1/2 

We recognize in (3.38) the Riemann-Hilbert problem of the orthogonal polynomials for the weig ht e-f dC- 
Specifically, if we denote by pI^\C,) the monic orthogonal polynomials that satisfy 



then the solution of the RHP (3.38), (3.39) is simply given by 

Rk = 5-^"^'iJK(C), 

where 



(3.42) 



(3.43) 



Hk{C) ■■= Co 



-2iTT 







VK-l 

-2z7rr^-2^ 



1 

2i7r 

-1 



P^k\^)^ 



* ds 



VK-l 
2iTT 



s 

P^j^\s)c-'""^" ds 



c 



Pi^Us)e 



C 



2^ + 2 

ds 



T:^{CoN)'' (3.44) 



VK-i " ^ ■ ■ VK-i Jr s-C 

It is crucial to point out here that the right multiplier TV"'*'^'^^ is needed to satisfy the correct jump relations (3.38), 
while the left multiplier Cq''^'^^ N'^^°'^ is needed to restore the boundary condition (3.39). One can satisfy the 
boundary condition (only) by choosing K as the closest integer of x. Defining S := x — K G ("ji 5)1 obtain 
the following estimate holding uniformly on the boundary. 



Rk 



0{N 



-7+2|<5|7 



), z e 



(3.45) 



Lastly, ^ kRk is analytic in D, because z~^'^^ factor in (3.43) cancels out the singularity of ^ k (4.2). 

The important observation is that \i x G 5 + ^ then the error term in (3.45) does not tend to zero (it is 0(1))- 
It is understandable as these values separate regimes where the value of K jumps by one unit and the whole strong 
asymptotic must changes its form. A similar problem arose in [3]. In section 4 we will overcome this obstacle. 

Remark 3.4 The orthogonal polynomials we are using here are a particular case of the so-called Freud Orthogonal 
Polynomials [16]. 



3.5 Asymptotic solution for Y and error term 

Collecting the results of the above analysis we have the following asymptotic solution for Y . 

k{z) outside of the disks 
l^as := { ^k{z)Rk z eD 

^'k(z)^ inside the disks around the turning points 



(3.46) 
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To find the error term we define the error matrix as follows. 



Y'i' j^^{z) outside of the disks 
YR^^'^K^iz) ZED 

YA^^'^](^{z) inside the disks around the turning points 



(3.47) 




The error matrix solves the 
residual RHP with the jump ma- 
trices as shown in the figure 8. 
It follows from the construction ^5>. 
that there is no jump inside the 
green disks, and on the cut. On 
the disks with airy parametrix 
and on dD, the jumps converge 
to the identity with uniform er- 
ror bounds in n L°° . We eval- 
uate these error bounds now. 

The jump of £ on is evaluated as follows. 



0{N- 



-7+2|(5|' 




V. 



Figure 8: The jumps of £. 



-7+27|5|-j 

The jumps of £ on the disks of Airy parametrix is similarly evaluated as follows. 



(3.48) 



(3.49) 



Therefore, the error matrix has jump matrices that are uniformly close to the identity with the error bound of 
C)(Ar-7+27l<5|) in L2 n A well-known theorem [15] guarantees that error matrix itself is bounded by the same 
error bound, i.e. £ = 1 + 0(Af-T+27|5|). This gives the following error term for the (strong) asymptotics of Y. 



Y = Y^sil + 0{N 



-7+27l<5| 



))■ 



(3.50) 



3.6 Necessity of improved approximation: a nonlinear Stokes phenomenon in x 

From the above estimate of the error term it appears that our global parametrix does a poor approximation if >c is 
not exactly an integer, and it is no approximation at all if >f G 5 + ^ since the error term is not vanishingly small. 
These transition points are the equivalent of the Stokes' lines in the standard theory of asymptotics of ODEs, where 
two solutions of the same RHP become of the same magnitude whereas off the line one is recessive and the other 
dominant. 

The phenomenon is similar here: for x e (^K — ^, K + ^) the dominant solution is the one we have constructed 
with K, where for x £ [K + + ^) it is the "next" with K +1. 

For X = X + 1 it is to be expected that both solution are of the same magnitude in a suitable sense and hence 
a "linear combination" should be sought for. 

We show how to do this in the next section and we will construct hence a parametrix including the subleading 
term in the N^"' expansion so as to have a uniform approximation to within 0{N~'^'^). 
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4 Asymptotic solution for x G M up to 0{N 
4.1 A short digression on Schlesinger transformations 

We start by observing that k and k+i are related by a Schlesinger transformation, as well as the local para- 
metrices Rk, Rk+i- In particular, we consider the "raising" Schlesinger transformation that raises the order of 
poles by one 

^k+i{z)=(^1+^^^k{z) , (4.1) 

for a z-independent matrix Sk- To show this, it is sufficient to observe that '^'k+i'^k^ has no jumps, behaves as the 
identity at infinity and has at most a simple pole at the outpost. The formula follows immediately from Liouville's 
theorem. Note also that Sk is a nilpotent matrix as follows immediately from the fact that det(l + Sk/ z) = 1. 
For future reference we compute Sk- Let us write the outer parametrix k as 

vI/;^ = [Aa',Bk]z^-% (4.2) 

where [A^, B^f ] is a 2 x 2 matrix holomorphic at z = 0. We also require det[Aif , B/f] = 1 to have det k ~ 1- 
Let us compute Sk in the raising Schlesinger transformation; the condition that determines Sk is that 

Ak(z) + -SkJ^k{z) = 0{z) (4.3) 
z 



or, equivalently. 



Ak(0) + SkA',,(0) = det[A,f(0),Ai.(0)l ' ' > ' 



Similarly the inverse "descending" transform 

^K = '^K-l, (4.5) 



Sk 
z 



requires the analyticity condition 



This determines Sk as 



^k{z)+^-Sk^k{z) = 0{z). (4.6) 
z 

~ det[V,BK(0)] , . , . 

^-^-det[B,(0),Bi,(0)]^-(°)' 

for an arbitrary vector V. It is also important to note that, since we have explicitly constructed the sequence of 
{^'/fjiir = 0, 1, • • •} in (3.19), the Schlesinger transformations between them must exist, and hence 

det[B^(0),B'^(0)] ^ ^ det[AA'(0), A'^(O)] , V/f e Z. (4.8) 

Now let us consider the similar transformation for the local parametrix. Previously in (3.43) we have constructed 
the local parametrix Rk = z^^^'^^Hk where 



K 



1 -uk/z 
-Ik-iIz 1 



l^^K (4.9) 
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All the components of the matrix in RHS have the multiplicative error bounds of 1 + 0{N '^'^ z ^). The two 
variables uk and tx-i are given by 

r^K TV^-r^-T 



UK 



2iTT Cq'^^+'i 



2^n Cl^^-- 

If- — 



(4.10) 



Note the interesting relation uk^k = 1, which will be essential for the consistency of our solution. 

From a similar argument used for deriving (4.1) using Liouville theorem, one obtains the following relation 



H 



K+l 



Z UK 

-£k 



H 



K- 



This transform can also be derived from the three-term recurrence relation 



CPKiO ^ PK+liO 



-PK-liO- 



The transformation matrix in (4.11) is LDU decomposed as follows 

L^'z'^'Uk, Lk := 



Z UK 
-iK 



1 

Ik/z 1 



Uk ■■= 



1 Uk I z 
1 



(4.11) 



(4.12) 



(4.13) 



From this we may view the matrices Lki Uk as the "two halves" of the transformation that raises the order of the 
Freud's OP by one. Therefore, these objects will appear for the "half-raising" transform, or more generally for the 
continuous transform that is parametrized by 6 in the next section. 



4.2 Improved parametrices 

Here we improve both the outer parametrix k and the local parametrix Rk to produce a legitimate asymptotics 
at half-integer x and to produce a better asymptotics for all x. Especially, we will use certain transformations 
that resemble the Schlesinger transformation discussed in the previous section. 

Looking back at the error analysis, the dominant error (3.45) originates from the off-diagonal term of Hkz^'^^ 
(4.9), especially from the terms u/z and 

A natural way to correct this problem is to define a new local parametrix,^ 



UkLk^iHk = 1 + 0{N-^''), zedD 



(4.14) 



so as to cancel out the leading off-diagonal terms. We get an improved error bound which is independent of x. 

We may also change the order of Lk-i and Uk in the above definition; although the two matrices Lk-i and 
Uk do not commute, the non-commutativity is within the error bound of 0{N^'^'') and, therefore, the order does 
not make any difference in the asymptotics. 

Given the local parametrix i?,^, we will find the corresponding improved outer parametrix ^'j^. It is notable 
that is constructed out of Lk-i and Uk, which appear in the raising transform (4.13). This suggests that 
the appropriate outer parametrix ^ ^ may come from a "partial 2-step Schlesinger" transformation of ^'k, such as 
^ = [1 + Fi/z + F2/z'^)'^K- Fi and F2 are then determined by the analyticity of ^R^ at z = or, equivalently, 

F-\ Fo 

+ — + 4 

Z Z'^ 



^kz~''^'UkLk^i=0{1), 



m 



(4.15) 



Clearly the issue is to remove the possible poles at z = 0, and the problem is addressed in the next section. 

®Our notation is not exactly consistent, since fj^ ^ R[f when >t = K. However, we believe that this will not make any confusion. 
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4.3 Improved outer parametrix 



The improved outer parametrix, that we will denote by ^'j^'', must satisfy the analyticity condition (4.15). is 
uniquely solved by the condition (4.15), and we will solve it in two steps by writing by first writing 



1 



K- 



Remembering k = [Aa', Jz^"'^ (4-2), we can rewrite ^R^ as follows. 



[Aa.Ba] 



1 uk/ z 
1 



1 

tK-l/z 1 



H 



K- 



Note that the terms u/ z and 1/ z depend on z and not on z. Accordingly we have the following definitions. 



[Ax,Ba] [Ak,Bk]{z/1^ 



Kl73 



z z 

1 



which is analytic at z = 0, and 



1 



K- 



(4.16) 



(4.17) 



(4.18) 



(4.19) 



The first step is to determine F by imposing the analyticity condition: (and we define another notation as 
below) 

1 Uk / z 



[Aa,Ba'] := 1 



F 



[Aa,B 







1 



= 0{l). 



(4.20) 



Assuming we solved the above, the second step is to determine G by imposing the new analyticity condition: 



[Aa,Ba] 



1 

^K-lIz 1 



0{l). 



(4.21) 



The above two analyticity conditions uniquely determine F and G, which can be written as follows using an 
arbitrary vector V. 



FV 



MAdet[y,AA(0)] 
l + UKdet[AK(0),A^(0)] 



Aa(0), GV: 



-^A-idet[V,BA(0)] 
l+^A--idet[B'^(0),BA(0)] 



Ba(0). 



(4.22) 



Here we have used the fact that det[AA(0), Ba(0)] = det[AA(0), Ba(0)] = 1. 

From the above formulae it is easily noticed that F ^ 0{uk) and G ~ ©(^a-i)- This immediately tells that 
[Ak,1^k\ — [Aa:,Ba-] + 0{uk) from the definition (4.20). This is useful to know because now we can change all 
the Bif's into B^-'s in the above equation (4.22) while keeping the error under O(ma^a-i) — 0{N~'^''). Since 
our asymptotic solution will have an error bound of 0{N~'^'') any term of that order or lower is meaningless. In 
addition, it is useful to observe the following facts for a further simplification. 



Aa(0) = Aa(0), det[AA(0),Ai^(0)] =det[AA(0),A'^(0)], 

Bx(0) = Bif(0), det[Bx(0),B'^(0)] =det[BA(0),B^(0)]+e)(MA). 



(4.23) 
(4.24) 



'^Note -once more- that ^ ^ ^ x when x = K 3.S for R^. However, we believe that this will not make any confusion since we are 
now constructing a refinement of the previous setting. 
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As a result we obtain (the first column of) up to 0{N as follows 

■■= (l + - + ^Kz'' (4.25) 



z I \ z 



~A , Mgdet[A;^,Ag(0)] cnw-ff-i det[A;^, Bk(0)] ^ mw-'f-i j. nc^^-27^ 

"^-^ +l + .,.det[A,^(0),A',(0)]^^^('^)" -l + £,_,det[Bi,(0),B^(0)]^-('^)^ + 

The second column is analogously given. In the last formula, when 5 ^ ±1/2, the first term provides the leading 
term, which becomes the strong asymptotics for OPs away from the outpost. Therefore, the strong asymptotics is 
5 independent. 

The subleading term is either from the second or the third term depending on the value of 5; to decide, we 
must recall the scaling behaviors of uk ~ TVt^^-'-i) and £k-i ~ iV-'^^^.s+i) (4.10). 

At both 5 — ±1/2, the leading strong asymptotics changes. At (5 = 1/2, uk is no longer scaling with N and, 
therefore, the second term also contributes to the leading asymptotic behavior. 

Mj I -a ..a- I MKdet[AK,Ag(0)] mw^-i 4^ror/V-27^ (Aoa\ 

*w|(i)-A;.z +l + „^det[A^(0),A'^(0)]^^(°)" ^^'^^^ 

One also observes that the number of roots at the outpost is still iiTasfor —1/2 < 5 < 1/2. 

At (5 = —1/2, Ik-i is the non-scaling parameter (see (4.10)) and, therefore, the third term contributes to the 
leading asymptotic behavior. 

\M I -A ^A-idet[Aj<-,BK(0)] ^K^ij, mw^-i 4-ncw-27^ (ao7-\ 

*A-i/2|(i)-A^z -i+^^_^det[B'^(0),BK(0)]" ^""^^^^ + ^^'^^^ 

In this case the number of roots is K — 1; one less than what it is for —1/2 < 5 < 1/2. The location of the missing 
root can also be found using the above expression. 
To summarize, 

i) At half-integer >c, the number of roots at the outpost is given by the closest integer that is smaller than 

ii) At half-integer >c, the outer parametrix ^ cannot be obtained by approaching from either side of >f. 

4.4 Roots at the outpost 

So far, we have described the strong asymptotics away from the outpost. Now we turn our attention to the inside 
of the disk D to look closely at the locations of the roots at the outpost. To this purpose, we evaluate ^'>ti?>c up 
to 0{N~'^'^) using (4.17). It is a straightforward but long calculation if one tries to obtain the full asymptotics 
up to 0{N~'^'^). Instead, here we will obtain only the leading term and the subleading term. (We write the full 
asymptotics in Appendix C for reference.) 

Looking at the first column of 'i'^R^ one gets 

(N-'C^f^^R^l^^ - (^Ak{0) + O{N-^+^^'\^))pP{0 (4.28) 

-(BKiO)+OiN-'>+^\'\'<))c^N'<eK^,pl,'UO- (4.29) 

Note that the above is the sum of two Freud's OPs. One may object that, say for a positive S, the second term is 
within the error of the first term, and cannot contribute as a subleading term. (For a negative S it is the first term 
that provides the subleading correction.) A closer look shows however that both error terms are C-independent up 
to 0{N~^) (which is not difficult to see from the general structure of the formula). So the error terms only change 
the coefficients of the two polynomials up to 0{N~'^). 
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From (4.28) we now identify the asymptotic locations of the 
roots at the outpost. Let us first consider a positive S. The roots 
are determined by the zeros of P^\0 the leading order. Ac- 
cording to the explicit value of [A/f,Bj<-] (3.25), the subleading 
term N~ 

-257pM^^((^) contributes with a same sign as the leading 

term. 

Due to the well-known interlacing property of the OPs, this 
means that all the roots are shifted from the zeros of P^\C,) by 
the amounts that scale as N^^^'' in (^-coordinate. 

In the figure 9 we show a schematic view of OPs and the roots 
where the real OP is shifted (to the left) from the leading asymp- 
totics (Freud's OP; dashed hne) by 0{N-^''^). 

For negative S the leading asymptotics is now given by Pi^liiC) 
and therefore we only see K — 1 roots to this order. Also from 
the exphcit value of [A^^ , Bj^], the subleading term N^'^"'^^^ P^\(^) 
contributes with the same sign as the leading term, which means 
that all the is: - 1 roots shift to the right by 0(7V257). 

Most interestingly there appears another root (let us call it "the 
stray root") to the left of K — 1 roots, distanced by ^ As 
in the schematic view (figure 10) the stray root scales differently 
from all the other roots. Though this root escapes to infinity in the 
(^-coordinate, it is actually converging to z = 0, only much more 
slowly than the other roots. It is also interesting to observe that, 
from the direction of the stray root, it seems to come from the 
main cut. The approximate location of the stray root is given by 





7\ 




n 












1 
















Figure 9: For S > 0, a, schematic plot of poly- 
nomial at the outpost. The dashed line is the 
plot of P^\(). The deviation between the 
two plots is of the order N~^^^ . 




]\1-2t\S\ 



Figure 10: For 6 < 0, a, schematic plot of poly- 
nomial at the outpost. The dashed line is the 
plot of P^'li(C)- One of the zeros is found 
away from the rest. In the z-plane, however, 
all the zeros converges to z = 0. 



Cstray 



(4.30) 



Finally let us consider the cases 6 — ±1/2. For both cases we get an additional subleading term. As the 
formulae are not particularly illuminating we present them in Appendix C. 



4.5 Kernel at the outpost 

Questions regarding universality hinge on the behavior of the kernel for the correlators [7] in the scaling coordinate 
In general the kernel is defined by 



1 Pn{z)pn^i{z') - Pn-l{z)pn{z') 
Z — z' 

1 [y-i(zOy(z)]2i ^ 

2j7r ^ ^ ~ 2:' 

J_ \y-^{z')Y{z)\-2,x ^_ N 

2i7r z — z' 



(4.31) 



19 



Using our asymptotics: Y ~ ^> ^R^, we evaluate the kernel in the local coordinate near the outpost 



(g(.)+g(.')+^) det «'^i?^(C)L.,«'^i?>.(C')L, 



7V27x(ArCo)^ r,^_i(C-C') 

Using the expression of ^'jti?,^ in the appendix we can in principle obtain the kernel up to 0{N^'^'') accuracy. Here 
we only show the leading term. The above approximation is valid for x ^ N + ^. This is exactly the kernel for 
Freud's OPs, as we expect from the previous sections. 

At K = K + 1/2 we obtain a different kernel at the leading order. 

"^"(^ ' ^) = - iV27><(ivco)7 m-i{(:-Q') ^ ^^^^ ^ "^^^ ] ^ ^ ^ 



where the constant ax is given by 



(4.34) 

:= ^^det[A^(0),A- (0)] 



1 + UK det[AKiO),A'j,{0)] 
At X = K — 1/2 we obtain the following kernel at the leading order. 

K(c'n e-f(^(-H^(-')+^) ( pP{OPi,'Un~p'K\C)p'KUc) , rnVi+t9rA.-^^ 

(4.36) 

where the constant f3x~i is given by 

^ ^K-idet[B^(0),BK(0)] . . 

^""^'^ l + ^K-idet[B'^(0),B^(0)]- ^ -^'^ 

4.6 Some consistency checks 

So far we have seen that the various asymptotic properties are drastically changed when x is a half integer. 
For X = K + 1/2 there are two ways to construct the asymptotics; one can apply a suitable (Schlesinger type) 
transformation starting from k or from In this section, the consistency of the two approaches will be 

proved in a completely general setting; without referring to the explicit global parametrix. 
When S = 1/2, we can obtain the same error bound using only Uk without Lk-i 

Rk+1/2 z-^'^^UkHk- (4.38) 
Correspondingly, the outer parametrix will then be written as 

«'k+i/2 := (4.39) 

where F has been explicitly solved for in (4.22). 

The same state can be approached from X + 1 by applying the "half-descending" transformation. In this spirit, 
the local parametrix may be written as follows. 

RT^,1:=z-''''+'^"^LkHk+,. (4.40) 
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Consistency means that -R^^^2 — Rk+i/2i which is an elementary consequence of the relations (4.11) together 
with the "factorization" (4.13). One realizes that uk^-k = 1 is the key identity. 

We can then find the corresponding outer parametrix ^'^^^2 demanding analyticity on '^^k+^/2^^k+i/2' 
prove the consistency of the outer parametrix we write 



K+l/2 



F 

z 



A'+l - 



F 

z 



Sk 

z 



(4.41) 



using the property of the Schlesinger transform. 

Therefore, the proof of consistency amounts to the following identity: 



F 

z 



To prove this, remember the following identity 



[Ak,^k\ = [AA'+i,B/f+i]; 



(4.42) 



(4.43) 



Here we do not use above A and B as we are dealing with a generic case. Then F and F are defined by the 
following conditions. 

1 uk/z 



1 



[Ak,B 







F 



[Ak+1, Bx+i] 



1 

1 

Ik/z 1 



= 0(1), 



= o{i). 



To prove (4.42) we only need to show that the following quantity is analytic. 



F 



Sk 



[Ak,Bk] 



1 uk/ z 
1 



0(1) 



1 

-Ik/z 1 



1 Uk / z 
1 



(4.44) 
(4.45) 

(4.46) 



It is straightforward to see that the above does not have a pole (hence analytic), using the identity uk(k = 1- 



This concludes that ^ 



(new) 



as it should. 



'K+l/2 - •^K+l/2, 

Now let us look at the kernel obtained in the previous section. At >f = A' + 1/2 the kernel has an additional 
term to the usual kernel from OPs. To make sense of the additional term we recall the following general identity 
for OPs. 



PniQPn-liC) ^ Pn{C)Pn-l{C) 



J=0 



(4.47) 



This tells us that the "raising operation" for the kernel is to add a term of the form oc Pn{C)Pn{C)i which happens 
to be the new term appearing in (4.34). 

As the kernel at x = K + 1/2 can be approached from either Kth kernel by adding aKPK{C)PK{0 o^' from 
K + 1th kernel by adding (3kPk{C)Pk{C)- The consistency of our result means the following identity. 

(P^k\OpP-.{C)~pP{C)p'kUO , aK 



vk{C-C') Vk 

'4^ii(c)4^Hc')-4^ii(c')4^^(0 , pk 



PKiOPKiO 



+ ^Pk{C)Pk{C) 



(4.48) 
(4.49) 
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This is reduced to showing aK — Pk ~ 1 which is easily obtained from the identity ukU-k — 1 and 

det[AK(0), A'^(O)] = — — -J— — (4.50) 

det[B'^_^^(0),Bif+i(0)] 

which follows from a direct computation using the Schlesinger transformation. 

There remains only one very exceptional case to consider: indeed both formula? in (4.22) may fail if the 
denominator appearing there vanishes. This can happen since at >f = if ± ^ one of the parameters u does not 
tend to zero. 

However such unlucky situation does not occur with our construction in (3.25). A direct computation shows 
that 

det[BWi(0),B;,^i(0)] = ^j^^^.i^V+i ■ (4-51) 



Plugging into the denominator of (4.22) we have 







l + ^^det[Bi^+i(0),Bif+i(0)] ^l + 2nT^ det[B'^+i(0), B;f+i(0)] (4.52) 

Vk 

Since rjK is the norm square of the monic Freud polynomial and given the sign of the imaginary part of (4.51), the 
denominator is strictly positive. 

4.7 Arbitrarily improved error bound 

Here we explain how to construct the outer and the local parametrix that has arbitrary small error bound as one 
wants to achieve. Though we will not explicitly carry out the evaluation of the corrections, the method already 
yields quite interesting identities which are otherwise hard to see. The main idea is to generalize the Schlesinger 
transformation to a higher order. 

Remark 4.1 The general framework for arbitrary improvement of the error is not new and appeared in ([17], Sec. 
7.2), based on the inversion of an operator close to the identity in terms of a Neumann (geometric) series. The 
approach of [5] to the problem was indeed based on those ideas. In a certain sense our approach is a manipula- 
tion whose "philosophical" meaning is the same as computing the terms of the above-mentioned Neumann series, 
although the practical details may be different. 

The Weyl function W{C,) is defined by 

^(0-^/1^--^-^- El- (4-53) 

(Let us write z instead of the correct z as we are not going to deal with the physical coordinate.) The measure was 
given by d/i(^) = exp(— ^^''+^) d^, but it can be general in the following discussion, {/ij} are the set of numbers 
defined by the expansion around z = oo. Note that they are also scaling with TV as /ij oc N'^^'y^^^ , 
Then the matrix, 

1 N^^fWiC) 
1 



(4.54) 
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satisfies the jump condition (3.38) for the local parametrix. Since the jump property remains by any multiplication 
to the left, we may multiply a matrix to get: 



1 Ej=iMi/2^' 

1 



1 iV2«7p^((^) 

1 



(4.55) 



which can be made as close to the identity (on the boundary of the disk) as one wishes by increasing R. Especially 
at i? = 2K we get the most modest error bound A^^t+^t^ ^ which we obtained in (3.45). 

The same error bound can be obtained in a different way, using the following property of the Weyl function 
and the Pade approximation provided by the orthogonal polynomials associated to the measure dfi{^) 



where 



1 

2j7r 



(4.56) 
(4.57) 



is a polynomial of order fc — 1 . 

Then the following matrix also has the same error bound as R2k- 



R 



K 



1 N^'^'^Qk{C)IPk{0 
1 



1 7V2«7p^((^) 

1 



-7+275 



(4.58) 



Now let us find the corresponding outer parametrix. Here it is necessary that we start from ^'q which does not 
have any pole at the outpost. Defining ^'o := [Aq, Bq] we propose the outer parametrix of the form: 



For ^nRn to be analytic near z = we demand 

[Ao,Bo] 




[Ao,Bc 



(4.59) 




1 Ej=iA^iA 
1 



Oil), 



(4.60) 



which completely determines {Fj}. As shown in the appendix, they are given by the solution of a linear equation. 

Now that we have explained the method to obtain an arbitrarily good error bound (on the boundary of the 
disk around the outpost) let us deduce a few implications. 

As we have seen already, the leading outer parametrix is given by the iV-independent "^k] viceversa we have 
just defined a set of TV-dependent outer parametrices ^'^ for R > 2K that should all converge to 'i'K, i.e. 




R > 2K. 



(4.61) 



Remember that {Fj} are determined by (4.60). 

Moreover, looking at (4.55), the local parametrix does not contribute to the first column of the full asymptotic 
solution "^fiRfi since Rk (or Rj^) are upper-triangular. Therefore, all the information about the asymptotics of 
the orthogonal polynomials near the outpost in this setting is encoded directly in the outer parametrix itself (4.59). 
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5 (New) Universality behavior 



Wc can now examine the result we have obtained in the neighborhood of the outpost. 

In particular we want to point out the behavior of zeroes of the the orthogonal polynomials Pn{z) in this 
particular scaling regime. 

• The normalized counting measure of the zeroes {z[^\ . . . , zi^^} of Pn{z) 



n 

^"(^)^=-E^(^-4"^) (5-1) 



converges in the sense of measures to the usual equilibrium measure, namely for any continuous function 
fix) e CO(M) 

/ u,,{x)f{x) = -j2 - / f{x)p{x) dx (5.2) 

• The fine behavior of those zeroes is precisely that in the scaling s {CoN)~''{x — ^o)- For < S < 1/2 
they are converging to the location of the zeroes {s[^\ . . . , s^'} of the orthogonal polynomial of degree K 
of the measure e~^^''*^ ds. In measure-theoretical terms, for any compactly supported continuous /(s) 

r ^ 

/ nUr,{x)f{{CoNy{x - Co)) = E/(4''')- (5.3) 

• For —1/2 < S < they are converging to the location of the zeroes {si^~^\ ■ ■ ■ ,s^]^Si^} of the orthogonal 
polynomial of degree K — 1 of the measure e"*^*^^^ ds. (There is also the unique "stray zero" that scales as 
^stray ^ ''^^^'^'''^^stray However, the constant Sgtray universal.) 

• the correlation functions are the same -in the same scaling- as for the random matrix model of size K 

dii{HK) := e-^^^K^'^ dHK (5.4) 
In a certain picturesque sense, there is a "microscopical" matrix model in the macroscopic background. 

• The kernel for the correlation functions Kn{x,x') = l^x^x']^^^^"'^^ ^ be computed from 

2nT X — x' 

and thus a direct computation gives the new universal kernel near the outpost (in the scaling coordinate): 

At S = ±1/2 there appears a new term - pI^\C)Pk\C) or P^}_i{C)PK-iiO depending on the sign - 
appears. 
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As we see, rather surprisingly, a finite-size matrix model (of particular type) arises naturally as a scaling limit of 
a general one. 

To conclude we remark that using the same methods employed here it is possible to handle a similar multicritical 
phenomenon where the "microscopic" matrix model has any polynomial exponential weight rather than the Freud 
one. This requires, however, a more finely-tuned potential, namely we need to introduce a dependence on A'^ into 
V and T. Similar considerations, but leading to Freud weights as in the main body of the text, are contained in 
Appendix A. 

We postpone this analysis to a future publication. 

A Double— scaling approach 

In this section we explain why our simplified approach using the chemical potential is de facto equivalent to a more 
refined double-scaling limit. In [6, 5] the deformation of the g-function was written in terms of the parameter T 
by fine-tuning T — Tcr {Tcr being the critical total mass of the measure): in [6] a small dependence on V was also 
introduced so that the critical point in the effective potential remains of the same order of vanishing, 2v + 2. Here 
we want to illustrate how our simplified approach can be related to those. 

Suppose that V :— V^{x),T :— depend on an external parameter which we denote by e in such a way that 
there exist a (finite) disjoint union of bounded intervals J = UJ^ with endpoints smoothly depending on e, a point 
^o(£) ^ J (also smoothly depending on e and a real function h{x) with the properties 

• h{x) is harmonic C \ J and continuous in C and at infinity h(x) ^ In \x\\ 

• ^V{x) - h{x) = for X £ J; 

• ^V{x) — h{x) = C{e){x — Co)^'^^^(l + 0{x — ^o)) — /(e), with /(e) a smooth function in a neighborhood of 
e = with M 3 /'(O) ^ 0; 

• other than the negative sign implied by the previous bullet-point in a neighborhood of (for small e) the 
sign of V{x) — h{x) on M \ J is strictly positive (see Fig. 11). 

• The sign of ^^V{z) — h{z) is negative on a left/right neighborhood of each component of J, and the size of 
this neighborhood is uniform in e for small e's. 

In the above C > may depend on e as long as it is smooth and bounded away from 0, and all the Lan- 
dau symbols should be uniform in e. The function h here is nothing but the real part of the ^-function (up 
to addition of the Robin constant) and the last bullet-point is equivalent to saying that h is the logarithmic 
energy of a positive measure supported on J (a consequence of Cauchy-Riemann equation for harmonic functions). 
In this framework then we could repeat verbatim the analysis by fine-tuning e via the implicit equation 

InA^ 

fie) = 2x7— (A.l) 

This equation defines e{N) (for suitably large N) since /(e) has nonzero derivative at e = and hence invertible 
near 0. 
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All the construction would still apply verbatim with 
the caveats that the turning points are (slowly) moving 
in N and hence the outer parametrix, the local coor- 
dinates ( used to define the local parametrices (includ- 
ing the one at the outpost ^q) depend (smoothly and 
slowly) on . 

The situation is not dissimilar to the one in [18] 
where an A^-dependent (^-function (with A'^-dependent 
turning points) was employed. None of the analysis we 
have carried out is significantly affected. In particu- 
lar the local RHP for the parametrix at the outpost is 
identical with the only understanding that the confor- 
mal parameter ({z) depends on A^ but only through e 
and hence in a uniformly bounded way. 

To convince the reader that the above list of require- 
ments is not insanely restrictive we show how to con- 
struct one such family of potentials V and total charges 
T. For simplicity we restrict to a one-cut situation but 
this is purely in the interest of conciseness. Methods similar to [13] could be used to construct a family for an 
arbitrary number of cuts. 

Suppose that V{x) is a (real) polynomial of even degree and that Vq^Tq is a multicritical pair as the one used 
in the main text. It is not difficult to show that we must have in general deg V >2v + A. 

We now define a deformation depending on the parameter e = T — T^r- Here we have chosen for transparency 
of exposition the deformation parameter as the deviation from the critical total charge, but in general it may be 
an abstractly introduced parameter. 

As explained for example in [13] we can write the (complex) effective potential ip^ as 




Figure 11: An example with v = \ and V of degree 6. 
This is a numerically correct plot, although the axis are 
scaled differently for V and ip. 



^o(z) = / Mq{x){x - S^of^+^^ix - aa){x -ba), < bo < 



bo 



where deg Mq = deg V' — 2^ — 2 and it is determined by the (algebraic) equation 

Mo{z){z - ^o)"'+V(^-ao)(^-M = K(^) - - + 0{l/z') 



(A.2) 



(A.3) 



For simplicity in what follows we assume that the roots /ij.o of Mq{z) are simple. It is clear that Vq cannot be a 
totally arbitrary potential (a simple parameter counting confirms this). 
The criticality condition amounts to 

r?o 

Mo{x){x - ^of^+Wi^ - '^o){x - &o) = (A.4) 

bo 

which implies that Mq{x) has an odd number of zeroes in (6, ^o)^- Other than this (unless other critical phenomena 
occur) Mo{x) > on [a,b] and Mo(^o) > 0. 



*By "smoothly and slowly" we mean that they are smooth functions of e which -in turn- is of order In N/N 
^Since degAfo > 1 the assertion on the minimal degree of V follows. 
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Define now 

Po{z) := Mo{xf{x - eo)'"+'(x - ao){x - bo) : (A.5) 

we will define a first-order ODE for with Cauchy data Pe\e=o = ^o- Since we want to preserve the multiplicity 
of the root at the outpost we must have — {z — + ...) and hence must be at least a polynomial 

of that order; moreover, since it should also preserve the square of M^, it must be divisible by AI^. The simplest 
choice is to take directly 

:-A-iM,(z)(z-C)^''+^ (A.6) 



The constant A is determined by the requirement that res = T = 1, namely 



^=°° 2^(z - a){z ~ b) 



A (A.7) 



This determines ^ as a (rather cumbersome) polynomial in ^. a, b: note that since we have chosen the deter- 
mination of the squareroot that is positive on the real positive axis near oo, then ^ < 0. 
The potential Ve undergoes the evolution according to 

V:{z)--+0{z-'). (A.8) 



2y/I\ ' Z 

Hence the coefficients of up to degree will necessarily depend on e. Only for the simplest case of = we 
can keep the potential fixed. 

The full ODE is thus either eq. (A.6) or -in terms of the position of the zeroes of P- 

" {A^+2)p^'-+'\o ~ W+mW^^W^^) ' 



P'{a) M{a){i-a){a~b) ' P'{b) M{b){i - b){b - a) 

PML (A 9) 

- 2P"(;.,) ^^-'^ 

The theorem of existence for ODE guarantees that the above equation has solution for e in a suitable interval 
around e = (with Pq as IVP). 

The only point we verify in addition is the claim about the behavior of the effective potential near the outpost 
^; but this is an elementary application of Taylor theorem since (noting that V'e(^e) = 0) 



^.i^) - / as+^^^^^h/KEMER^, _ + oix - 0) (A.10) 

The derivative of /(e) at e = is 

m = ^r d.<o (A.11) 

^0 Jbo v(a;-ao)(x-6o) 

This implies that for e > (T > Tcr = Pq) the effective potential is negative in a small neighborhood of the 
outpost ^, with a minimum value that should be fine-tuned as detailed above. 
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B Construction of the outer parametrix for arbitrary number of cuts 



We only sketch the construction since the details would require a good deal of notation to be set up. We will use 
the same notation and ideas contained in [12]. 

We denote by w double-cover of the z-plane branched at the endpoints of the support of the equilibrium measure 

-.^ Uiz-a,) (B.l) 

This is a hyperelliptic algebraic curve of genus g. We denote by oo± the two points above z = oo in the usual 
compactification of the curve, and by p± the two points projecting to the location of the outpost^". We denote by 
ujj the first-kind differentials normalized along the a~cycles: explicitly 

^ dz 

ujj{z) = aji — - — , (B.2) 

where the summation over repeated indices is understood (and they range from 1 to g) and aj^ is an invertible 
matrix such that ujj = djk- 

Using the standard notation for divisors on Riemann surfaces [19] we consider the unique (up to multiplicative 
constant) sequence of spinors with the divisor properties 

{i;i°^)>-{r-l)oo+ + roo- + K{p+-p_), reZ (B.3) 
4°>{p):^4'Hp*) (B.4) 

where p p* is the holomorphic involution of the hyperelliptic curve. The spinors (and their starred counterparts) 
are also sections of the Une bundles £,£^^ with character x (x~^ respectively) defined by 
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e""'^-^^ ■.= Y[e'^'' for7 = aj 
e=i 



2 /"""^ 
f 



/ p{x) dx = j-th filling fraction (B.5) 



'Q2j-1 

These are the generalization to arbitrary genus of the spinorial Baker- Akhiezer vector used earlier: the char- 
acteristics A,B gC^ are defined up to integers. 
The matrix 



(B.6) 



dz 

solves the model RHP with quasi-permutation monodromies on the cuts and on the gaps (antidiagonal on the cuts, 
diagonal on the cuts). Here p^^{z) is the point {z,W-^-{z)) where w+(z) is determination of w{z) that behaves like 
z^^^ at infinity, analytically extended to the complex plane sliced along the support of the equilibrium measure 
(i.e. the physical sheet). 

The spinor V dz is defined on the double cover of the hyperelliptic curve (it has \/ branch-behavior at the 
Weierstrass points in terms of the local parameter — aj , hence has singularities of type z — a when thought 
of as a spinor on the plane) . 



^"The point oo+ is characterized by lii > as 2 S M-|_ near oo. The point is the point on the Riemann surface of w obtained by 
analytic continuation of w on the complex plane slit along (a2j — i,02j) 



28 



The first column has a zero of order K 
at — Co (the outpost) and the sec- 

ond column a pole of order (at most) K: at 
infinity it behaves as z'^'^^ (up to left multi- 
plicative constants). 

The expression in terms of Q functions 

is 



Vr,K '■ = 




"29+2 





A 



(p+rcx)_ - {r+l)co+ + K{p+ - p-)) /^^(p) 


Qr+l 


oo+)ef (p-p_) yfdz 



Vr,K{p) 



K 



Gr K = lim 



Vr,K{p) 

C'r,K = lim 

p^oo_ z^ 



£ a.^ 6,(0)0.,- 



(B.7) 
(B.8) 
(B.9) 

(B.IO) 



The notation (rather standard) is lifted from [12] and [20]: the Abel map is understood when writing points as 
arguments of Q and it is based at one of the Weierstrass points (for example ai) A is an arbitrary odd non-singular 
half-period. Recall that (pag. 23 of [20]) all such characteristics A are in one-to-one correspondence with partitions 
of the Weierstrass points into g — I and g + 3 points {aki , . . . , afeg_i} U {o^^j ' • ■ ■ > "^fe It is to be noted that (pag. 
23 of [20]) 

hA{p) 



v/nLi^(2-"fc* 



and that 



e 



(p + roo_ - (r + l)oo+ + K{p+ - p_)) /ia(p) 



F{z) 



Q^{p- oo+)Vdz 



n -If (--«.) 



(B.ll) 



(B.12) 



where F{z) is an analytic function with jump discontinuities on the cuts^^ and it is independent of the choice of 
A, it is bounded in the finite region of the z-plane and growing like z^^ at infinity. Moreover, straightforward 
computations show that (with some overall ambiguity of signs) 



e,(p-oo±) - T-die^{0)aeg 

p-*oc± z 



hA{p) 



CK,r — 



^ , ^ ~- (Ta,e,(o)a,,)^5 

Q^(p- oo±)V dz p^oo± 

e;(oo+ - oo_)ef (oo+ - p+)e \^] (r(^_ - + k{p+ -p_)) 



Ck.t — 



(a,e,a,<,)'^+ief (oo_ - p+)e f^l ((r + 1)((^_ - 00+) + k{p+ -p_)) 



e^+i(oo_ -oo+)ef (oo_ -p_) 

Of course this is a bit vague description since we should first stipulate how the fourth-roots have been defined. 



(B.13) 
(B.14) 

(B.15) 

(B.16) 
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The constant Cr,K is simply the normahzation so that V'r behaves as z"^ at oo-i- while Cr^K is chosen so that 
ipr{p*) behaves as z ^ ^ at oo _i_ . 

Remark B.l The fact that these formal expression do not vanish identically follows from the fact that 



e 



(r(oo_-c5o+)+i^(p+-p_))^0 (B.17) 



for arbitrary A (z , as these correspond to positive divisors of degree g with g points in the gaps. This is proved 
in a more general setting in [Prop. 6.3, pag. Ill of [20]]. Of course for K = = r the nonvanishing of this 
expression is precisely the same that appears in [2] although maybe not clearly stated. 

From a point of view of isomonodromic theory, the above theta function is intimately related to the isomon- 
odromic tau function [21], whose vanishing determines the (non) solvability of a Riemann-Hilbert problem. In this 
example the RHP is the model problem for the OPs. If the spectral curve had no real- structure then in general it 
could happen that for exceptional values the problem does not admit a solution (see [12]) 

We conclude this section with a few important remarks and shortcomings of these formulas 

• the construction of the improved parametrix in Sect. 4 did not use the specific form of the outer parametrix 
but just the jet-expansion near the outpost, thus it applies verbatim to the general case, with the proviso of 
the next point; 

• the description of the behavior of the roots at the outpost remains generically valid in this case: however the 
direction of approach of the stray zero (in eq. (4.30)) depends on the actual sign of the expressions involved, 
hence in this general case it cannot be easily identified^'^; 

• for the Stokes' values >f G Z + ^ there is the potential for the denominators of formulae (4.22) to vanish under 
exceptional circumstances (i.e. for special spectral curves and special values of K). This would make the 
approximation (4.25) unbounded in N and hence invalidate it. For the one-cut case as in the main text it was 
rather simple to directly verify that the determinants in (4.22) have a suitable sign so that the denominators 
are bounded away from zero, but for the case of multi-cut solutions a similar computation requires a deep 
manipulation of functions and we could not determine a similar property. We suspect that such property 
should hold here too on account of the reality conditions of the cuts and the Jacobian of the spectral curve. 

C Asymptotics: long results 

We write the full asymptotics for the first column of Y up to 0{N~'^'^). 

{N^C2f^^R^\^,^ (C 

= P^"VO^A I UKdet[AK,AKm Ak(0) ^^^i det[Ag,B;^(0)] 6^(0) \ 

\ " l + UKdet[A,^(0),A'^(0)] ^ 1 + £^-1 det[B'^(0), Bk(0)] ^ ^ ) 



Nor it should be expected to always come from one side. On a heuristic level, the stray zero should come from the "closest" 
spectral band, and hence it depends on the location of the outpost. 
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W ^ / V l + £K-idet[B'^(0),BK(0)] ^ y 

-^x-iiV^C^P*'^^ (0 (bk+ ^^det[BK,B^(0)] Ak(0) _ det[BK,Bj^(0)] 6^(0) ^ 

° V l + WKdet[AK(0),A'^(0)] z 1 + ^^.i det[B'^(0),BK(0)] z z ^ 

+0(7V-2t). 

At half integer x we have the following leading and the subleading behavior. 



l + ^x;fdet[A;^(0),A'^(0)] V 

(iV Goj M/.K.|(,,_ i + ^^_^det[B'^(0),BK(Or^-^^^^ 

R' mi / det[A'^(0),B;,(0)]-[z/i]i tn.\pH(.s 
(^A^(O) + Ik-MO) - iK-i i + ^^_^dct[Bi,(0),BK(0)]^^^°) j 

/ ^RWm ^K-idet[B'^(0),BK(0)] \ (.) 

-^K-i ^B,(0) - i + ^^_^det[B',(0),B,(0)]^-(°V ""^^^^^ ^ 

for X = iV' + 1/2 and x = -fC — 1/2, respectively. 
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